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Chiral Dynamics and S-wave contributions in Semileptonic D s /B s decays into tt + ti 

Yu-Ji Shi 1 *, and Wei Wang 1,2 1 
1 INPAC, Shanghai Key Laboratory for Particle Physics and Cosmology, 

Department of Physics and Astronomy, Shanghai Jiao-Tong University, Shanghai, 2002f0, China 
2 State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 

Chinese Academy of Sciences, Beijing 100190, China 

In this work, we study the semileptonic decay modes B° —> iv + n~£ + £~ and D+ —> 7r + 7r -£+ v 
in the kinematics region where the 7r + 7r _ system has a invariant mass in the range 0.5-1.3 GeV. 

These processes are valuable towards the determination of S-wave 7r + 7r _ light-cone distribution 
amplitudes whose normalizations are scalar form factors. We compare the results for scalar form 
factors predicted in unitarized chiral perturbation theory and extracted from the data on the B a -A 
J/%j)7 r + 7r - . Then the B s —> n + n~ and D s -A n + n~ form factors are calculated in light-cone 
sum rules, based on which predictions for differential decay widths are made. The results are in 
good agreement with the experimental data on the B s and D a decays into 7r + 7r - . More accurate 
measurements at BEPC, LHC and KEKB in future will be helpful to examine our formalism and 
constrain the input parameters more precisely. 


I. INTRODUCTION 


It is anticipated that new physics (NP) beyond the standard model (SM) can be indirectly probed through the 
precision exploration of low-energy processes. An ideal platform is to study the flavor-changing neutral current 
(FCNC). Rare B decays like the b —> s £ + £~, with tiny decay probabilities in the SM, are sensitive to NP degrees of 
freedom and thus can be exploited as indirect searches of these unknown effects. In terms of observables ranging from 
the decay probabilities, forward-backward asymmetries, polarizations to a full angular analysis, the exclusive decay 
mode B -A K*£+£~ can provide us with a wealth of information on weak interactions. Recent measurements of the 
almost form-factor independent ratio P 5 by the LHCb collaboration have indicated a deviation from the SM by about 
3.7a 0,3. 

In fact, the B -A K*£ + £~ is a four-body decay process since the K* meson is reconstructed in the Kn final state. 
Thus it is more appropriate to explore the B -A M\M 2 £ + £~ , in which various partial-waves of Mi M 2 contribute 0-Q|. 
The S-wave contributions to B —> Ktt£ + £~ have been discussed for instance in Refs. 0 Q- The bottom mass mb is 
much heavier than the hadronic scale Aqcd> which allows an expansion of the hard-scattering kernels in terms of the 
strong coupling constant a s and the power-scaling parameter Aqcd /mb- On the other side, final state interactions 
among the two-light hadrons should be constrained by unitarity and analyticity. A formalism that makes use of 


these two advantages has been developed in Refs. 12Ml4j. and summarized in Ref. 15J. Such an approach was 


pioneered in Ref. flbl . [r7 ] , and a method without the analysis of hard-scattering kernels has been explored recently in 
Refs. 00- See also Refs. 0-0] for attempts to analyze charmless three-body B decays. The aim of this work is 
to further examine this formalism by confronting the theoretical results with the relevant data on B s -A /i - 

and D s -A 7 r + 7 r - e + v e . 

In Ref. 0|, the LHCb collaboration has performed an analysis of rare B s decays into the 7 r + 7 r - p + p~ final state 
with the measured branching ratio: 


B{B S -a n+TT-fi+pL-) = ( 8.6 ± 1.5 ± 0.7 ± 0.7) x 10" 


(1) 
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where the errors are statistical, systematic and arise from the normalization, respectively. The dominant contribution 
is the B s —> f 0 (980)fj, + fj,~ [ 28 ]: 


B(B S / o (980)(—)■ tt+tt" )/x+/x") = (8.3 ± 1.7) x 10 


1-8 


29 


( 2 ) 

3- 


This has triggered theoretical interpretations based on two-meson light-cone distribution amplitudes (LCDAs) 

Previously, the CLEO collaboration has investigated the D a —> n + n~£u e , in which the /o(980) contribution is found 
dominant as well 


31 


321 


B{D S / o (980)(-l> Tr+Tr")e+i/ e ) = (2.0 ± 0.3 ± 0.1) x 10 


-3 


A recent analysis based on the CLEO-c data |33j] gives a similar result 

B(D S —s- / o (980)(—>■ TT+Tr“)e+i/ e ) = (1.3 ± 0.2 ± 0.1) x 10" 3 . 


(3) 


(4) 


The BES-III collaboration will collect about 2/6" 1 data in e+e" collision at the energy around 4.17GeV, which will 
be used to study semileptonic and nonleptonic D s decays [341. 

The rest of this paper is organized as follows. In Sec. El we will give the results for scalar tttt form factors and the 
non-local LCDAs. Sec. Ell will be devoted to the calculation of the B s —> tt + tt and D s — > tt+tt form factors in the 
light-cone sum rules (LCSR). In Sec. IIVI phenomenological results for a variety of observables in the B s —> tt + tt~£ + £~ , 
B s — > tt+t x~vv and D s — ^ tt+tt ~ £v are presented, and compared to the experimental data if available. An agreement 
between theory and data will be shown in this section. Our conclusions will be given in Sec. [V] 


II. SCALAR tt+tt" FORM FACTORS AND S-WAVE LCDAS 


A. Scalar form factor 


We start with the definition of a scalar form factor: 


(tt+tt |ss|0) = B 0 F^(m nn ), 


and the Bq is the QCD condensate parameter: 


< 0 |ot| 0 > = -/=fl 0) 


(5) 


( 6 ) 


with /tt as the leading order (LO) pion decay constant. For the numerics, we use f n = 91.4MeV and (0|qq|0) = 
— [(0.24± 0.01)GeV] 3 [35j |. This corresponds to B 0 = (1.7 ± 0.2) GeV. 

In the literature, a variety of approaches have been used to calculate the ^^(m^), including the (unitarized) 
chiral perturbation theory (yPT) 36-0 and dispersion relations 44|. In the xPT, the LO and next-to-leading order 
(NLO) results can describe the low-energy data with a good accuracy, but with the increase of the invariant mass, 
higher-order contributions become more important. This is under expectation since the perturbative expansion in 
XPT is organized in terms of p„/(4 tt f n ). It has been argued that the unitarized approach, by summing higher order 


corrections, can extend the xPT applicability to the scale around 1 GeV [45|]. A fit to the BES data on the tttt 
invariant mass distributions in J/ip — > tt+tt~<)> S has been performed in this approach and an overall agreement is 
found 0| . The fitted result for F^(rn^) is shown in Fig. [U The modulus, real part and imaginary part are shown 
as solid, dashed and dotted curves, from which one can observe the broad structure for the cr(600) and the peak at 
/o(980) is naturally produced. 

In recent years the LHCb collaboration has conducted a series of analyses of angular distributions in the B s —> 
J/t/ttt+tt" decay mode 47, 48'] . In this procedure, the S-wave contributions have been explicitly separated, and three 
resonances, /o(980), /o(1500) and /o(1790), have been identified. To access the tt+tt" invariant mass distribution, 
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FIG. 1: The in unitarized xPT. The modulus, real part and imaginary part are shown in solid, dashed and dotted 

curves. The result is based on the fit of the BES data on the 7T7 t invariant mass distributions in J/tf) —> Tv + Tr~(j> 0] in Ref. [3tj |. 


the Breit-Wigner formula is employed for the /o(1500) and /o(1790). Due to the fact that the /o(980) lies in the 




vicinity of the KK threshold, the Flatte model [43, [50( has been adopted. Considering the relative strengths and 
strong phases among different resonances, we have 


KA™ 2 ™) = 


c i m / o (98O) ei01 


rn n 7T m /o(980) ^ m /o(980) 


9kkF]^ k pkk) 


C 2 m / o ( 1500) e * e2 


~ m / o (1500) + * m /o(150O) r / o (i5O0)(^7i-) 

_ C3TO /o(1790) et6>3 _ 

m Tnr — m /o(1790) * m /o(1790)F/o(1790) (to^*.) 


The p*.*. and pkk are phase space factors 


47 


48 


50(: 


Ptttt - „ \ / 1 


4m!j 



1 - 


4m; 


Pkk — 



(7) 


( 8 ) 


Compared to the normal Flatte distribution, an additional correction has been introduced in the LHCb fit above the 
KK threshold to better describe the data 


48] 


Fkk = exp (-ak 2 ), 

where k is the kaon momentum in the KK rest frame, and a is set to a = 2.0GeV~ 2 
width r s(m^) for an .S'-wave resonance is parameterized as 


(9) 


48| . The energy-dependent 


r sKJ = r s 


ms 


TO* 


— 4m; 


m| — 4?7l^ 


to 


( 10 ) 


with the constant width Tg, and the Blatt-Weisskopf barrier factor Fr = 1 [48(. The c* and 6i , i = 1, 2, and 3, are 
tunable parameters. In the fit by LHCb, two solutions are found and the fitted parameters for three contributing 
components are collected in Tab. HI [48]. 

In Fig. [2[ we compare the results for the scalar form factor calculated in the unitarized %PT (dashed curve), and 
the ones extracted from the B s data based on Eq. © (solid and dotted curves). The constant c\ in 

Eq. © has been tuned in the comparison. The two solutions found by the LHCb give very similar shapes as shown 
in the figure. There are a few remarks on the shapes. 
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TABLE I: Fitted parameters for contributing components in the B s — > J/npir+it by the LHCb collaboration [48|. Two solutions 
are found in the fit. 


Fractions (%) 

Solution I 

Solution II 

/o(980) 

70.3 ± 1.5+5 t 

92.4 ±2.0^5 * 

/o(1500) 

10.1 ±0.8to.3 

9.1 ±0.9 ±0.3 

/o(1790) 

2.4 ± 0.4l* ° 

0.9 ± 0.3±o j 

Phase differences (°) 

Solution I 

Solution II 

/o(1500) - / 0 (980) 

138 ±4 

177 ±6 

/o(1790) - /o(980) 

78 ±9 

95 ±16 

Parameter 

Solution I 

Solution II 

mf Q ( 9so) (MeV) 

945.4 ±2.2 

949.9 ±2.1 

g-nn (MeV) 

167 ±7 

167 ±8 

Q K K / Q'K'k 

3.47 ±0.12 

3.05 ±0.13 

^n/oCisoo) (MeV) 

1460.9 ±2.9 

1465.9 ±3.1 

P/ 0 (i 500 ) (MeV) 

124 ±7 

115 ±7 

m/ 0 (i790) (MeV) 

1814 ± 18 

1809 ± 22 

P/o(i790) (MeV) 

328 ± 34 

263 ± 30 



FIG. 2: The predicted in the xPT (dashed curve), and fitted from the B a — » J/ipn+ir data (solid line for solution 

and dotted line for solution II). The two fitted solutions are not distinguishable below 1.2 GeV. 


On the contrary with the yPT results, the parametrization in Eq. 0 does not contain the / o (500) (or the 
so-called cr) contribution. The LHCb collaboration has set an upper limit for the ratio 0|: 


F r (fo(500)) 

Wo (980) 


< 0.3% 


( 11 ) 


at a 90% C.L. 


• The expansion parameter in yPT is p v / (4n f^), where the p n is the pion’s momentum. Summing higher order 
s-channel contributions and incorporating the coupled channel effects will extend the applicability region up to 
lGeV. On the other hand, the parametrization in Eq. 0 has explicitly included two scalar resonances, /o(1500) 
and /o(1790), and thus is valid in the region above 1 GeV. 


At the /o(980), the unitarized yPT leads to a very narrow peak. This feature is smeared out in the data since 
the binned results for the 7 T 7 t invariant mass distributions are presented on the experimental side. Still the 


parametrization in Eq. 0 is well consistent with the data on the B s —> J/ipTr + n 


47, 


481. 
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• From the comparison, we can see that at the current stage the advantages in both cases are not overwhelming. 
In future, wc hope the results can be improved in some more sophisticated methods like the unitarized approach 
with resonances [Tf . In the following, we will use the parametrization form inspired by the data in Eq. 0. 


B. Generalized LCDAs 


With the scalar form factor as the normalization condition, the S-wave tt+tt LCDAs are defined by 


12 


51 


54|: 


((tt+tt ) 5 |s(j:)7 m s(0)|0) = F° n (m nn )p WWtlJ , / due mp ^' x <j) nn (u), 

Jo 

((7r+7r-) s |s(a:)s(0)|0) = F^(ml n )B 0 [ due iup ™' x ^(u), 

Jo 

1 f 1 

(('K + n~)ss(x)a^s(0)\0) = - ywa^) / due mp ^' x (j)^{u). (12) 

O Jo 

The LCD A qb nn is twist-2, and the other two are twist-3. Their normalisations are given as 

f ducf) s nn (u) = [ du^(u) = 1. (13) 

Jo Jo 

The conformal symmetry in QCD [ 55 )] indicates that the twist-3 LCDA have the asymptotic form 51 -Q: 

&n( u ) = 1 > CttM =6u(1-u), (14) 

while the twist-2 LCDA can be expanded in terms of the Gegenbauer moments: 

•(«) = 6u(l — u) a n C^/ 2 (2m — 1). (15) 

n 

In most cases, contributions from higher Gegenbauer moments are suppressed and thus one may keep the lowest 
moment ai. It is worthwhile to stress that these LCDAs for a two-hadron system have the same form as the ones for 
a light scalar qq meson 51 -0. In Ref. 0, the first Gegenbauer moment for the /o(980) is calculated as 1 


ai = -1.35. 

while the perturbative QCD analysis of the B s decays has used a much smaller value (29f : 

ai = —0.36. 


(16) 


(17) 


III. HEAVY-TO-LIGHT FORM FACTORS 


The B s -A tt+tt transition can be parametrized by three form factors 


((tt+tt )s\s'y^ 5 b\B s ) = 


((tt+7t )s\sa liv q ,/ 'y 5 b\B s ) = 


—1 

m Bs 

T l 

J rr 


P»- 


m r - m: 


-9m 


-r-B s —>-7nr / 2 2\ , 

Pi {m nw ,q) + 


m r - mz 


r T > . 


B s —>mr /■ 2 




VL,? 2 ) 




(m| s - ml^qn - g 2 P^ 


(18) 


1 In Ref. [5(J, the normalization factor for a scalar qq meson is m / 0 (980)// 0 (980) > which is the F^^Bq in this work. The results for the 
twist-2 Gegenbauer moment in Eqs. m and Eli have been converted to our convention. 
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where the orbital angular momentum in the 7r + 7r~ system is chosen as zero in order to select the S-wave. The p Bs and 
p n7r is the momentum for the B s and the tttt system, respectively. The momentum transfer is defined as q = p B — P 
and the P is defined as P ^ = p Ba + Pnn- Here the convention slightly differs with the ones adopted in Refs. 

The D s —> 7r + 7r _ form factors can be analogously defined, with the replacement ma, —> m Bs - 
As we have demonstrated in Ref. 


li- 


7 T+7T LCDAs 


51 




1J], the LCSR allows us to express the B s —> tt + tt form factors in terms of the 


5J, |57H59| . T le LCSR factorization formulas read as [121 ]. 


= n f 


du 


-exp 


m 2 + uumA — uq 2 


+ 


ml + q 2 - u 2 ml n 


uM 2 


uM 2 
+ exp 


-^Mu) + uK n (u) + ^(u) 


Jo 
' M 2 J 


KMo) m l - u l m A + q 2 

6 ml + - 9 2 


(19) 


PI 


\ml^q 2 ) = 


I f ^ 

Nf s / —exp 

Uuo u 


mt 




uum^ n — uq 


uM 2 


m b 

~Bo 


$ ™(u) + (2 - u)<f>l n (u) 


u{ml + q 2 - u 2 mA) + 2 (m 2 b - q 2 + u 2 ro 2 J ®° n ( u ) 


3 u u 2 M 2 

Uo(jn 2 b + q 2 - u 2 mU + 2 (m\ - g 2 + 
uoiml + ulm 2 ^ - q 2 ) 


6 


exp 


s o 

' M 2 J 


^>o) 

6 


( 20 ) 


Pf^im^q 2 ) = P, 


B s —>n 7T 


(ml n ,q 2 ) + 


-P: 




( m W7T, 9 2 ), 


( 21 ) 


P%“^ n (mA,q 2 ) = 2N F (m Ba +m nn )\ [ —exp 


(ml — uq 2 + uum 2 n ) 
uM 2 


+m b 


^%(up) exp [-sp/M 2 ] 

6 ml - q 2 + ulmh, 


®nn(u) , „ KnA) 
2B 0 b 6uM 2 


( 22 ) 


where 


Nf = B 0 FA(mlA Tl + T‘ exp 


2 m B jB a 


M 2 


uo = 


min + q 2 -SQ + \J (m 2 ^ + q 2 - s 0 ) 2 + 4 ml n (mj -g 2 ) 
2m 2 


(23) 


In order to derive the above equations, the Borel transformation of hadronic and of QCD expressions of correlation 
functions has been carried out, defined as: 


B[P(Q 2 )} = lim 


Q 2 — »oo, n— >oo, (yi ^I 


1 HQ 2 r(-L) HQ 2 ) 


dQ 2 


(24) 


where J 7 is a function of Q 2 = —q 2 and M 2 is the Borel parameter. The explicit form is: 


B 


1 


exp(— s/M 2 ) 


(s + Q 2 ) n \ (M 2 ) n (n - 1)! ' 


(25) 


This operation improves the convergence of the OPE series by factorials of n and, for suitably chosen values of M 2 , 
enhances the contribution of the low-lying states to the correlation function. 

For convenience we can define the normalized form factor: 


Pi(mA,q 2 ) = B 0 FlA m in)Fi(q 


( 26 ) 
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FIG. 3: At the maximal recoil q 2 = 0, the dependence of F\{q 2 = 0) = Fo(q 2 = 0) (left panel) and Fr(q 2 = 0) (right panel) 
on the Borel parameter M 2 . The final results are obtained requiring stability against variations of M 2 . 


TABLE II: Fitted parameters of the B a /D a —» 7r + 7r form factors derived by LCSR. 


B s 7T + 7r 

Fi(q 2 = 

0) ai 

bi 

D a n + 7r 

Fi(q 2 = 

0) ai 

bi 

Fi 

3.66 

1.39 

0.54 

F i 

2.45 

0.82 

0.20 

F 0 

3.66 

0.54 

-0.08 

Fo 

2.45 

0.39 

-0.15 

Ft 

4.29 

1.33 

0.54 






where the r and q 2 dependence has been factorized into the (rn 2 ^) and Tk(g 2 ), respectively. This approximation 


is justified by the Watson-Madigal theorem |60|,|61). As a reference point, we will choose the m n7V = mf 0 ( ggo) to explore 
the functions Fi{q 2 ). 


In the numerical analysis, we use [62l 1631 

f B , = (224 ± 5)MeV, s 0 = (34 ± 2)GeV 2 , B 0 = (1.7 ± 0.2)GeV. 


(27) 


With these numerical inputs, the sum rules (fTHll - d^l) provide us with the functions Fi(q 2 ) for each value of q 2 as 
a function of the Borel parameter. In Fig. [3] at q 2 = 0 we show the dependence on M 2 with a± = —0.6. The 
results are obtained requiring stability against variations of M 2 . As demonstrated in this figure, the form factors 
become stable when M 2 > 12GeV 2 . The situations with different ai and q 2 values are similar and thus we can choose 
M 2 = (16 ± 2)GeV 2 . 

Since the form factors are sensitive to the Gegenbauer moment a± of the twist-2 LCDA, in the left panel of Fig. [4] 
we show this dependence in the range ai = (—1.4, —0.4) at the maximal recoil q 2 = 0. We will show later the value 
ai = —0.6 can describe well the data on both the B s —> ir + ir~£ + £~ and D s —> ir + Tr~lv. 

The LCSR is applicable in the hard-scattering region. To access the momentum distribution in the full kinematics 
region, we will adopt the following parametrization 


F t (q 2 ) = 


Fi( 0) 


1 -a l q 2 /m 2 B +b l (q 2 /m 2 B ) 2 ' 


(28) 


where i = 1,0, T. The parameters can be fitted in the region q 2 < 5GeV 2 for the B s transition, and the results are 
collected in Tab. GH 

For the D s —> transition, we use [63] 


so = (6.5 ± l)GeV 2 , f Da = (257.5 ± 4.6)MeV, m c = 1.4GeV. 


(29) 


















FIG. 4: At the maximal recoil q 2 = 0, the dependence of F^o , and Ft on the Gegenbauer moment a\ is shown in the left 
panel. Dashed and solid curves correspond to Fi(q 2 — 0) = Fo(q 2 = 0), and Fr(q 2 = 0), respectively. In the right panel, the 
q 2 dependence is given with a i = —0.6. Solid, dotted and dashed lines denote the Fr(q 2 ), Fi(q 2 ) and Fo(q 2 ), respectively. 




FIG. 5: The functions Fi(0) for the D s 


the dependence on M 2 (Gegenbauer moment ai) in the left (right) panel. 


The results for the D s —> tt + tt~ form factors are given in Fig. [5] From the left panel, we can see the results are 
stable when M 2 > 6GeV 2 , and we will use M 2 = (8 ± l)GeV 2 . The dependence on the first Gegenbauer moment 
ai is less severe compared to the B s —> 7r + 7r _ case, as shown in the right panel of Fig. [5] In the D s mode, the 
twist-2 contributions in the regions with x > 1/2 and x <1/2 cancel with each other. This fact has been explored 
in the study of D s —> /o (980) transition [64]. Since the energy release in the D s transition is small, we have used the 
—5GeV 2 < q 2 < 0 to fit the g 2 -dependent parameters in Eq. 







9 


IV. PHENOMENOLOGICAL RESULTS 


A. B s — > 7T + 7T V 


We proceed with the analysis of B s 
nian ,6a] 


r”t *rr 


£ + £ , whose decay amplitude is governed by the effective Hamilto- 


H eS = -%V tb V t * s Y, CiitiOib). 


10 


V 2 


i=1 


The Oi is a four-quark or a magnetic-moment operator, and the C)(/x) is its Wilson coefficient. The explicit forms can 


be found in Ref. 65J. Gf is the Fermi constant, and V tb = 0.99914±0.00005 and V ts = —0.0405tg;oi2 I63j] are the CKM 
matrix elements. The bottom and strange quark masses are m b = (4.66±0.03)GeV and m s = (0.095±0.005)GeV [631. 

In general, various partial-waves of two-hadron M 1 M 2 state contribute to a generic B —> £~ process and 

the differential decay width has been derived using the helicity amplitude in Refs. 0-Q. In the B 
the S-wave contribution dominate with the angular distribution: 


d 3 T(B s —> 7T + 7T /x + /r ) 

dm 2 7 T dq 2 d cos Qi 


J'i + J 2 C cos(20<) 


where the coefficients are 


7T+7T fl + fj, , 


(30) 


J°i = |^o| 2 + l^ol 2 + 8^I^Vflolcos(^ 0 -^ 0 )+4m,V? 


J°2 = -Ph Ko| 2 + l^ 


flOl 


(31) 

(32) 


In the above equations, P 21 = -y/l — 4m|, rhi = mi/\fq 2 1 and 61 is the polar angle between the B s and the /i moving 
direction in the lepton pair rest-frame. The d° 0 and <5^ 0 are phases of the helicity amplitudes 


A° , 

•^L/R, 0 




= N i 
~^°L/R,t = l ^N% l i 


1 


m Be 

1 

™-B s 


(C 9 T C w ) 


Vx_ 

V 7 ? 




■B s —tmr / 2 


0 q ) + 2(CVl — Cm) 


V\m b 


\fq 2 {m B + m^) 


Tr 


B s —>7nv (2 


(<? 2 ) 


(C 9 tC w )- 


m n — m: 


X = A% t - A° L . t = 2 N 2 \/n^C 10 1 




1 

m Bs 


7T7T -pB s —>-7T7r ( 2 
^0 


( 9 2 ) 


171 Bs -£ Bs —>7T7r ^2 ^ 




. 7T 
•J r 


0 


where 


N? = %—VM 


4v^ 7T 

1 


Nf = 

2 16tt 2 


V 1 -4ml/ml n x - 


8 V\q 2 /) 2 i 


3 2567r 3 m^ 


The Kallen function A is related to the 7T + 7r momentum in the B s rest-frame: 

A = A(mg s , m 2 +n -, q 2 ), A (a, b , c) = a 2 + b 2 + c 2 — 2 (ab + bc + ca) . 


(33) 

(34) 

(35) 

(36) 

(37) 


In Fig. [6l results for differential branching fractions dB/dm^n for the B s —> qi + fi~ are given in panel (a) 

and (b), and the ones for B s -A are given in panel (c). The result in panel (a) clearly shows the peak 

corresponding to the /o(980). In order to compare with the experimental data 0, we also give the binned results 
in panel (b) in Fig. [G] from 0.5 GeV to 1.3 GeV with square markers. Dominant theoretical errors arise from the 
B 0 = (1.7 ± 0.2)GeV. The experimental data (with triangle markers) has been normalized to the central value 
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FIG. 6: Differential branching ratios dB/dm^n for the B s —> n + n~fi + fi~ in panel (a) and (b), and B a —> n + n~r + T^ in panel 
(c). In panel (b), experimental data (with triangle markers) has been normalized to the central value of the branching fraction: 
B(B a — > n + n~ fi + ) = (8.6 ± 1.5 ± 0.7 ± 0.7) x 10 -8 0]i an( f theoretical results are shown with square markers. 



FIG. 7: The differential branching ratio dB/dq 2 for the B a — > n + n p (solid curve) and B s — > n + n t + t (dashed curve) 
is given in unit of 10~ 8 /GeV 2 . 


in Eq. ©■ The comparison in this panel shows an overall agreement between our theoretical predictions and the 
experimental data. Integrating out the n from 0.5 GeV to 1.3 GeV, we have the branching fraction: 

B{B S TT+Tr’/i+p") = (6.9 ± 1.6) x 10“ 8 , (38) 

which is also consistent with the data in Eq. ©■ 

In Fig. 0 we predict the differential distribution dB/dq 2 (in unit of 10~ 8 /GeV 2 ) for the B s — > n + n~/j, + /r _ 
(solid curve) and for the B s —> (dashed curve). Results for the integrated branching fractions of 

B s —> tt + -k~t + t~ are predicted as: 

B{B S —> 7T+7T -T+T~) = (8.8 ± 2.1) x 10“ 9 , (39) 


where 0.5GeV < m 
facilities including the LHCb detector 


< 1.3 GeV is assumed. Our theoretical results could be examined at the future experimental 
and the Super-B factory at the KEK 671. 
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»W(GeV) (a) 



FIG. 8: The differential branching ratios for the B s —► 7r + 7r vw. the left panel for the dB/dm nn , and the right one for the 
dB/dq 2 . 


B. 


B s 


The b —$■ svv effective Hamiltonian is given by 

G p O^em 


H b 


y/2 27rsin 2 (0w) 


V tb V* Vx X(x t )0 L = C l O l , 


which involves the four-fermion operator 


0 L = [s7 ft (l - 75)&][i / 7/i(l “ 75 W ■ 


(40) 


(41) 


Here 9w is the Weinberg angle; the function X(xt) {xt = m 2 /rn 2 v , with rrit the top quark mass and mw the W mass) 
has been computed in Refs. 65, 6 SJ, and the QCD factor ijx is found close to one 69- 7ll]. 


65|, | 68 j], and the QCD factor ijx is found close to one 
With the above Hamiltonian, we obtain the differential decay width 

d 2 T(B s — > ( tt + tt~)vv) 


= 3 x |Aq| , 


dq 2 dml n 

where the factor 3 arises from three species of neutrinos. The helicity amplitude in this case is 

1 7A 

1 


(42) 


HO = Cl^N^i— 




T^(m\^q 2 ) 


N 2 " = 

167T 2 


y /1 -4x 


8 V\q 2 


3 2567r 3 m^ s 


(43) 


We give our predictions for the differential distributions for the B s —> n + n~vv in Fig. [ 8 ] the left panel for the 
dB/dm^n, and the right one for the dB/dq 2 . The integrated branching fraction in the range 0.5GeV < m n7r < 1.3GeV 
is predicted as: 

B(B S tt + tt~vv) = (4.9 ± 1.2) x 10" 7 . (44) 


There is a large chance measure this branching ratio at the Super-B factory at KEK 67j. 


C. D s —> 7T + 7T tv 

The effective Hamiltonian for c —> stv transition is given as 

= N([s 7 m (1 - 7 5 )c][i> 7 M (l - 75 )t] + h.c., 


(45) 
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with 


The differential decay width for D s —> tt+tt~ £vt can be expressed as 

d 3 T 3 r 


dm 2 K7T dq 2 d cos 0/ 
with the li having the form: 


h(q ,m K7V ) + I 2 {q , m Kn ) cos(20 ; ) + J 6 cos (0 t ) 


(46) 


(47) 


h(q 2 ,ml n ) = [(1 + rhf)\A° 0 \ 2 + 2m?|A°| 2 ] , 

h{q 2 ,rnl n ) = -( 3 i\A° 0 \ 2 , 

J 6 (g 2 ,mL) = 4rh 2 Re[A 0 0 A° t *}. 


(48) 


Using the D s —> tt+tt form factors, the matrix element for D s decays into the S-wave rnr final state is given as 


= N‘JNZi- 


1 


where 


As discussed in Ref. [l 
dependent ratio: 


A'i = Nf, 


m = 


'Nit 


m DB 

1 

m Ds 


Va 




D s —>mr 


(m. 


2 Q 2 ) 

7T7T 5 1 ) 


m D — m 




7V7T -pD s -A-ITU 
^0 


(ml^,q 2 ) 


167T 2 


\/l - 4 ml/min x 


8 \/\q 2 


3 2567t 3 to|j 


(49) 


(50) 


one can explore a number of the g 2 -dependent ratios and in particular the lepton flavor 


7 Z+ /e {ml " 2 


<n = 


d 2 T(D s — > tt+tt /i + u ll )/dq 2 dm 


4* 


and the integrated form over g 2 : 


R^(m 2 J = 


d 2 T(D s —> Tr+TT-e+i' e )/dq 2 dm 2 7 


dr(D s —> TT+TT [i+v^/dr 


dT(D s —> 7r+7r e+v e )/dm. 


2 

7T7T 


(51) 


(52) 


Our results are given in Fig. 0 The first panel corresponds to the dB/dm^n, hi which the dotted and solid curves 
denote to the D s —> tt+tt~ fiv and D s —> TT+Tr~ev, respectively. One different behavior in the differential branching 
ratio with B s —> tt+tt~ £+l~ in the m W7r distributions is the suppression in the large m region. The second panel 
shows the ratio R^^ e (rn^. v ) defined in Eq. (15^1) . A comparison with the experimental data on the differential branching 
, 32| is given in panel (c), where we can also find the agreement. 


fraction 


31 


The integrated branching fractions are predicted as 


B(D S 

B(D S 


>TT+Tr~e+v) = (1.52 ± 0.36) x 10' 3 , 
tt+tt~h+v) = (1.68 ± 0.39) x 10“ 3 , 


(53) 

(54) 


where 0.5GeV< ro„ < 1.3GeV has been adopted in the integration. Again the errors come from the QCD condensate 
parameter Bq. Theoretical results are in good agreement with the CLEO results in Eqs. (13141) [3 lit! .'ll. We expect 
experimental errors will be greatly reduced since in future the BES-III collaboration will collect about 2/& -1 data in 
e + e _ collision at the energy around 4.17GeV which will be used to study semileptonic and nonleptonic D s decays [341. 
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FIG. 9: The differential branching ratios for the D s — > -r + tt~£v. The first panel corresponds to the dB/dm^n, in which the 
dotted and solid curves correspond to D s —> TT + n~ fj,u and D s — > 7r + 7r _ ei' respectively. The second panel shows the ratio 
defined in Eq. ((52)). A comparison with the experimental data [3l|, [32| is given in panel (c). 


V. CONCLUSIONS 

Rare B decays have played an important role in testing the SM, and hunting for the NP. In recent years, a lot of 
experimental progresses have been made on the B —> K*£ + £~, and remarkably the LHCb collaboration has found 
a 3.7 a deviation from the SM for the ratio Pg. This observable Pg is believed almost independent on the hadronic 
uncertainties. 

The analysis of B —> Vl + £~, more appropriately B —► M\M 2 £ + l~, requests not only the knowledge on the 
Wf , expansion but also the M\M 2 final state interactions. In this work, we have studied the P° — > tt + tt~£ + £~, 
B£ —> and Df —> t: + 'k~£ + v decay in the kinematics region where the 7r + 7r - system has a invariant mass in 

the range 0.5-1.3 GeV. These processes are dominated by the S-wave contributions and thus they are valuable towards 
the determination of the S-wave 7 t + 7t _ light-cone distribution amplitudes which are normalized to scalar form factors. 
We have compared the results for scalar form factors calculated in unitarized yPT and the ones extracted from the 
data on the B s —> J/ipn + n~. We have derived the B s —> and D s —> 7r + 7r“ transition form factor using the 

light-cone sum rules, and then presented our results for differential decay width which agree well with experimental 
data. Accurate measurements by the BES-III at the BEPC, the LHCb at the LHC and Super-B factory at KERB in 
future will be valuable to more precisely examine our formalism, and determine the two-hadron LCDA. 
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